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Abstract: Water resources are critical for the existence and development of oases in endorheic basins. Thus, 
to enable sustainable development, it is fundamentally important to understand how to allocate and use 
these resources in a reasonable way. We therefore simulated and analyzed changes in water consumption 
pattern within the Dunhuang Oasis of China under three scenarios using a system dynamic model that 
corresponds to different water consumption pattern. This was done to assess the impacts of regional water 
resource planning (comprehensive planning of the rational use of water resource and protection of 
ecosystem services in the Dunhuang Basin) on water consumption pattern within the Dunhuang Oasis. The 
first of these, Scenario 1, is a baseline in which the status quo is maintained, while Scenario 2 incorporates 
the comprehensive effects of agricultural water-saving irrigation measures with an inter-basin water 
diversion project, and Scenario 3 focuses on ecological rehabilitation. In the baseline Scenario 1, the total 
water consumption within the Dunhuang Oasis increased progressively while agricultural water 
consumption remained extremely high and threatened overall ecological security. In contrast, Scenario 2 
would decrease agricultural water consumption by almost 5.30107 m° following the implementation of 
water-saving practices. The additional water allocated from an inter-basin water diversion project would play 
an important role in alleviating ecological strain on the oasis. Finally, in Scenario 3, the total irrigated land 
must be decreased to 20.6107 hm? by 2025 assuming that water supply for ecosystem restoration would be 
at least 50% of the total consumption. Although water resource planning plays a very important role in 
alleviating the ecological water crisis within the oasis, it is necessary to consider the suitable scale of oasis 
with regard to current water consumption pattern. 
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1 Introduction 


Arid and semi-arid regions cover more than 30% of the Earth's land surface (Okin et al., 2006) and 
22% of land area in China (Hu et al., 2010). An oasis is defined as a specific landscape that exists 
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within deserts in arid and semi-arid regions and that enables flourishing vegetation and human 
settlement because of a stable water supply (Ling et al., 2013). Without water, an oasis cannot 
survive and will become desertified; this coupled relationship between an oasis and water is 
determined by basic water resource patterns. For a catchment in arid region, the upper reaches 
is the water production area with snow and glacier melting from mountains; and the middle and 
lower reaches are the water consumption area with scarce precipitation. Enhanced water 
consumption in the upper and middle reaches caused by the booming population and rapid 
economic development results in more serious ecological problems in the corresponding lower 
reaches (Cheng et al., 2014). These effects include the desiccation of terminal lakes, increased 
sandstorms and the death of Populus euphratica forests leading to a series of severe ecological 
disasters (Department of Geoscience, Chinese Academy of Sciences, 1996). The world heritage 
sites of Mogao Grottoes, Mingsha Mountain (sounding dunes), and Crescent Spring all are 
located within the Dunhuang Oasis of China (Zhang et al., 2014). The water resource crisis in 
this region is therefore a threat to this natural oasis as well as the precious cultural heritage. 

The rapid development of intensive agriculture alongside an increase in human activity, 
especially production expansion and the modification of irrigation channels in the middle 
reaches, have substantially modified the catchment hydrological cycle and have negative 
impacts on ecosystem sustainability within the lower reaches (Qi and Luo, 2006; de Fraiture et 
al., 2010). Therefore, the main problem is the competition in water demands between the 
developments of socio-economy and sustainable ecosystem within inland river basins (Micklin, 
1988; Feng et al., 2001). A huge range of water resource planning strategies have therefore been 
implemented to address this issue within inland river basins around the world, including the 
Tarim River Basin (Sun et al., 2011), Heihe River Basin (Cheng et al., 2014; Hu et al., 2015a), 
and Shiyang River Basin in China (Zhu and Li, 2014), the Murray-Darling River Basin (MDB) 
in Australia (Wheeler et al., 2014) and the Guadalquivir River Basin in Spain (Berbel et al., 
2015). 

The MDB Plan is an example of an early water resource planning strategy that was introduced 
in 1914 to develop the River Murray Waters Agreement following several years of droughts, 
debates and community actions to address water shortages (Crase et al., 2004; Woodhouse et 
al., 2006; Burke, 2012). The MDB Plan was finally signed into law in December 2012 (Zhao et 
al., 2015) and has since had a marked effect on recovering ecosystem functions and balancing 
the multiple needs within this basin (Bark et al., 2013). The most advanced water resource 
planning strategy currently in place in Europe is the case of the Guadalquivir River in Spain; 
this waterway is the largest river in southern Spain and crosses 25% of the nation's irrigated 
area. The Guadalquivir River Basin has reached its extraction limit, however, and has no 
additional water sources available in this region. A consequent reduction in water use is 
therefore planned across this basin, including water-saving measures to be applied in both 
agricultural and urban sectors. This plan represents a useful strategy to achieve sustainable water 
management because it details economic and technical measures to control demand while 
fulfilling wastewater treatment objectives including non-point source pollution from agriculture 
as well as urban point source pollution (Gutiérrez-Martin et al., 2017). The first successful water 
resource planning case within China has been the Heihe River Basin Plan (Cheng et al., 2014). 
The Heihe River Basin is the second largest inland river basin within the northwestern China. 
Excessive water consumption due to agriculture production has impacted artificial oases in the 
middle reaches, and the Ejin Oasis in the lower reaches has experienced increasingly serious 
ecological problems. An ecological water diversion project was therefore implemented at the 
end of the 20" century and severe ecosystem deterioration in the lower reaches has been 
dramatically mitigated (Cheng et al., 2014). The middle reaches of this basin are required to 
transfer 9.50x 10° m? of water to the lower reaches in "normal" years (defined as a year in which 
the upper reaches of the basin discharge 15.80x10°® m? of water) (Hu et al., 2015b). This 
increased water volume is then allowed to flow into the lower reaches of the basin restoring 
lakes and supporting the recovery of vegetation within the Ejin Oasis (Hu et al., 2015a). 

The focus of this analysis is the Dunhuang Oasis, which is adjacent to the Heihe River Basin. 
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This oasis has experienced competitive water pressures between economy and ecosystem, which is 
similar to that of Heihe River Basin. The Danghe River has a total volume of 4.63x10* m? and is 
the largest watercourse within the Dunhuang Oasis. The water resource exploitation and utilization 
rates are close to 100% in the oasis and the agricultural water consumption accounts for nearly 90% 
of total resource available. Ecological water use over recent decades has been impacted by the rapid 
development of the local economy. This has led to an ecological crisis mainly manifest in a 
shrunken terminal lake, a decline in the groundwater level and accelerated desertification. Water 
has become both the primary restrictive factor and bottleneck for socioeconomic development 
within the Dunhuang Oasis. Thus, in order to protect precious local heritage and improve the 
ecological environment, Dunhuang City has proposed a program of "Comprehensive Planning of 
the Rational Use of Water Resource and Protection of Ecosystem Services in the Dunhuang Basin" 
(DBP) (Water Resource Department of Gansu Province, 2011). Programs proposed within the DBP 
include agricultural water-saving measures, patterns of crop and industrial structural adjustments, 
as well as a water transfer project from the Sugan Lake Basin to the Danghe River Basin. The 2020 
goals of the DBP are to: (1) maintain groundwater in an equilibrium state; (2) raise the water level 
of Crescent Lake to 2 m to recover its original shape; (3) stop the expansion of desertification 
within the basin and grassland degeneration around the oasis fringes; and (4) maintain the West 
Lake Nature Reserve in a stable state. Scenario 2 and Scenario 3 in this study are based on the DBP. 
We also extended the simulation period assessed here from 2020 to 2025 to encompass the extra 
time necessary to implement this plan. 

This study therefore emphasized the DBP that aims to simulate future water consumption 
patterns under three scenarios. Baseline Scenario | is based on the situation at present, Scenario 
2 is formulated according to the DBP framework, and Scenario 3 incorporates ecological 
recovery. Analysis and comparison of these three scenarios provide a partial basis for policy 
making across this area as well as other dryland regions. This is important because the ecological 
problem of Dunhuang Oasis is not isolated but is common in other inland river basins as well as 
arid regions globally (Niu et al., 2010; Zhang et al., 2014). 


2 Study area 


The Dunhuang Oasis (40°02'-40°28'N, 94°31'-94°55’E) is bounded by West Lake Wetland to the 
west and Danghe River Reservoir to the south (Fig. 1). The elevation of this region is between 80 
and 1650 m a.s.l. As this area is far inland, a dry climate with low precipitation is typical. Based on 
data collected from 1954 to 2014 by Dunhuang Meteorological Station, the annual mean temperature 
within this area is 9.6°C with recorded maximum and minimum temperatures of 43.6°C in July and 
—28.6°C in January, respectively. Although the mean annual precipitation is only39.2 mm, the mean 
annual potential evaporation can reach 2486.0 mm. Precipitation is concentrated between April and 
September in this area, the average relative humidity is less than 40% and the annual solar radiation 
is 6418.4 MJ/m’. 
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Fig. 1 Location of the study area derived by Landsat TM image (Bands 5-4-2) 
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3 Methods 


The system dynamics (SD) approach was first proposed by J.W. Forrester who is known as the 
founder of SD. This method describes system structures using stocks, flows and internal feedback 
loop functions in diagrams to simulate dynamic behavior (Forrester, 1961). A SD model can 
therefore be used to describe linkages among system components at different levels and can also 
illustrate dynamic feedback processes. The SD model can also be used to predict changes in 
complex systems under "what-if" scenarios and can therefore be very useful for examining and 
recommending policy decisions in management and social contexts (Mohapatra et al., 1994). The 
SD methodology has been applied across numerous research fields including urban and eco- 
environmental planning as well as for the evaluation of land carrying capacity (Higgins et al., 1997; 
Gertseva et al., 2004; Huang et al., 2009). In this study, we established, tested and simulated water 
consumption pattern (WCP) under three scenarios using the SD model. 


3.1 SD model components 


On the basis of the relationship between variables and feedback structure, we determined a causality 
diagram and a flow chart (Fig. 2) and then used the software Vensim PLE SD (Ventana System, 
Inc. USA) to design a model. Under WCP, water supply comprises surface water, groundwater 
exploitation and upcoming inter-basin water diversions while consumption includes agricultural, 
domestic, industrial, service and ecological components. We simplified this model into four 
subsystems comprising water resource, population-livestock, economic and ecological 
components. The "what-if" scenario assumption applied here as well as the causal relationships 
between individual subsystems are described in more detail in subsequent sections. 
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Fig. 2 Components of the water consumption pattern (WCP) in Dunhuang Oasis 


3.1.1 Water resource subsystem 


Surface and groundwater are the most important components that comprise the water resource 
subsystem within the study area (Fig. 3). Indeed, when the Sugan Lake Basin water diversion 
project is completed, the Dunhuang Oasis will receive an additional volume of water about 
8.35107 m°. To simplify this model, water supplied from this project was assumed to be zero prior 
to 2020 and to then remain at a constant volume of 8.35x107 m? between 2021 and 2025. The 
diversion project aims to recover the ecological conditions of both the Crescent Spring and West 
Lake Wetland and to ease the severe groundwater depletion of Dunhuang Oasis. 
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Fig.3 Water resource subsystem in Dunhuang Oasis 


3.1.2 Population-livestock subsystem 
The primary domestic water consumer in this case is the population-livestock subsystem. This 
subsystem comprises urban domestic, rural domestic and livestock water consumption (Fig. 4). 
Specific variables and their relationships to one another are listed in Tables 1—6. The feedback of 
the population-livestock subsystem is reflected by the annual increase rate of water demand and 
the shortage of water available (Fig. 5). 
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Fig. 4 Population-livestock subsystem in Dunhuang Oasis. "+" that means the latter increases as the former 
increases or the latter decreases as the former decreases, which is defined as positive feedback. "—" means that the 
latter increases as the former decreases or the latter decreases as the former increases, which is defined as negative 
feedback. 


3.1.3 Economic subsystem 

The economic subsystem in this case comprises industrial, service industrial and agricultural water 

consumption. Industrial and service industrial water consumption volumes are influenced by 

corresponding industrial outputs, technological processes and water use efficiencies (Fig. 6). 
Agricultural water consumption is influenced by the changes in total irrigated area, crop patterns 

and water-saving practices (Fig. 7). Thus, considering that the DBP is feasible, this study assumes 
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Fig. 5 Population-livestock subsystem loop 


that: (1) irrigated area will decrease gradually by 3600 hm? until 2025; (2) high ecological value 
cash crops (i.e., grapes, fruit trees) will be developed while the number of food crops (i.e., spring 
wheat, corn) will be decreased; (3) irrigation efficiency will be improved by increasing micro- 
irrigation, canal micro-irrigation and greenhouse micro-irrigation areas to ensure that the entire 
Dunhuang Oasis becomes a water-saving region. These measurements aim to decrease net 
agricultural water consumption, replenish the growing domestic-industrial water demand and 
ultimately alleviate groundwater depletion. 
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3.1.4 Ecological subsystem 


Ecological water demand comprises urban and rural ecological water demand as well as ecological 
water discharge and the Crescent Spring ecological water supplement (Fig. 8). Urban ecological 
water demand includes the water requirements of urban public green and lake areas as well as urban 
roads. Rural ecological water demand, in contrast, primarily consists of the wind break forest 
system water requirement, while ecological withdrawal denotes the amount of water released 
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directly from the Danghe River Reservoir. The Crescent Spring is a world famous natural landscape 
that also requires a certain amount of water to maintain a stable water level. 
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Fig.7 Agricultural subsystem in Dunhuang Oasis 
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Fig. 8 Ecological subsystem in Dunhuang Oasis 
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3.2 Data 


The data used in this study include observation records between 1990 and 2010 obtained from the 
Dunhuang Meteorological Station as well as population and industrial records extracted from the 
Dunhuang National Ecological Statistics Yearbook (Bureau of Statistics of the City of Dunhuang, 
2007, 2013), water use quota data extracted from industrial water use records for Gansu Province 
(Water Resources Department of Gansu Province, 2005), from the field water-saving 
transformation project preliminary design report of the Danghe River irrigation district (Gansu 
Province Water Conservancy and Hydropower Investigation and Design Institute, 2011), and water 
resource planning data extracted from comprehensive plans for the rational use and protection of 
ecosystem services within the Dunhuang Basin (Water Resource Department of Gansu Province, 
2011). These data were used to calculate the service industrial, industrial and agricultural water 
consumptions as well as the ecological water demand. The statistical data used in this study did not 
include that of Yanguan Town. 

The specific equations and variables included within the model utilized here are summarized in 
Table 1; this summary presents the partial logic relationship between domestic and agricultural 
water consumption, as all other subsystems exhibit a similar relationship. Records in Tables 2—4 
include initial values and table function values and "if then else" functions. 


Table 1 Logical relationships among variables 


System item Variable Code Logical relationship 
Domestic Total population 11 11=INTEG(12,11) 
water Annual increasing amount of population 12 12=11xT1x(a+bx17/c) 
consumption Population growth rate T1 - 

Urban population 13 13=11xT1 

Urbanization rate T2 - 

Rural population 14 14=11-13 

Rural residents living water-use coefficient F1 - 

Urban residents living water-use coefficient F2 - 

Rural domestic water consumption 15 15=14xF1 

Urban domestic water consumption 16 16=13xF2 

Water shortage degree 17 - 
Agricultural Irrigated area 21 - 
water Cotton area 22 22=21xT3 
consumption Cotton area proportion T3 - 

Cotton canal irrigation area 23 23=22xT4 

Cotton micro-irrigation area 24 24=22xT5 

Cotton canal irrigation area proportion T4 - 

Cotton micro irrigation area proportion T5 - 

Cotton canal irrigation norm F3 - 

Cotton micro irrigation norm F4 - 

Cotton water requirement 25 25=23xF3+24xF4 

Note: T denotes a variable that is a table function; F marks a variable that is a "if then else" function. INTEG, integral; -, no logical 


relationship with other variables: a=1; b=2; c=100. 


Table 2 Table function values 


Variable 1990-2020 1990 1995 2000 2005 2010 2015 2020 
Urbanization rate —0.006—1.000 0.195 0.255 0.436 0.487 0.510 0.592 0.649 
Service-industrial water used 

(x10? CNY) 9.000—25.000 25.000 22.000 20.000 18.000 14.500 11.300 9.300 


Cotton canal irrigation area 0.000-1.000 1.000 1.000 1.000 0.970 0.880 0.710 0.610 
proportion 
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Table 3 Initial values 


Variable Initial value 
Total population (x10*) 10.99 
Large livestock inventory (x10*) 2.38 
Small livestock inventory (x10*) 48.08 
Industrial output (x10’ CNY) 3.28 
Service industrial output (x10’ CNY) 6.04 
Irrigated area (x10* hm?) 2.56 


Table 4 Summary of "if then else" functions 


Variable Table function 


If then else (Time<2001, 0.043, if then else (Time>2001: and: Time<2009, 
0.05, if then else (Time>2009, 0.065, 0))) 

If then else (Time<2000, 0.095, if then else (Time>2001: and: Time<2009, 
0.106, if then else (Time>2009, 0.13, 0))) 

If then else (Time<2009, 700, if then else (Time>2009: and: Time<2015, 
650, if then else (Time=>2015, 600, 0))) 

Note: These functions illustrate values during different time periods. 
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3.3 Model calibration and validation 


This study encompasses a simulation length between 1990 and 2025 with records at annual 
intervals. The total simulation period considered here spans two stages, calibration and prediction. 
The first of these, the calibration stage, spans the period between 1990 and 2010, while the latter 
spans the period between 2011 and 2025, simulating future WCP under different scenarios. 

Model calibration was conducted prior to analysis. In this step, simulation results were compared 
with actual historical data to verify the extent of their agreement as well as to assess the reliability 
of model parameters and the accuracy of the simulation model. Simulation results for total 
population (Fig. 9a) as well as for industrial and service industrial production (Figs. 9b and c) were 
compared with actual data. Records in Table 5 summarize comparisons between simulated and 
historical data for the period between 2006 and 2010; relative errors were considered acceptable 
when less than 5% absolute value (Table 5). Considering the randomness and variability of both 
the natural environment and social policy, these results suggest that our model is both valid and an 
effective simulation. 
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Fig.9 Simulated results for total population (a), industrial production (b) and service industrial production (c) 
3.4 Scenario simulation 


Once the model had been calibrated, WCP scenario analysis was performed. Scenarios were 
constructed to both forecast and simulate potential changes in watershed management. The 
development of the Dunhuang Oasis depends on subsystem coordination, including population, 
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environment, economic and ecological variables. In order to coordinate contradictions between 
resources, environmental and regional development provide the premise for the sustainable 
development of the Dunhuang Oasis. According to the local situations and actual needs, this paper 
simulates the WCP based on three scenarios, present-day as the baseline (Scenario 1), the water- 
saving/water-diversion condition (Scenario 2), and the ecological rehabilitation scenario (Scenario 


3). 
Table 5 Accuracy assessment of the simulated results 


Relative error (%) 


= Total population Industrial production Service industrial production 
2006 0.00 —0.01 1.09 
2007 0.00 —0.01 0.54 
2008 0.00 —0.01 4.04 
2009 —0.11 —4.71 3.75 
2010 0.38 —1.38 4.97 


Note: Yangguan Town was excluded from the statistical data presented. —, simulated value<real value. 


3.4.1 Scenario 1 


This scenario, also known as the baseline, assumes that development policies and system structure 
do not exert a large adjustment on the forecasting period. The records in Table 6 include all 
functions as well as "if then else" operations within the system. Thus, considering current 
development states, regional restrictive factors and policy implications, this research assumes a 
moderate trend function in terms of parameters. Urbanization rate was set at 59.0% for 2015 and 
65.0% for 2025, while the industrial and service industrial water consumption decrease from 160.00 
to 52.70 and 18.00 to 9.30 m? per 104 CNY, respectively, via the implementation of water-saving 
measures. Numerous functions exist within the agricultural subsystem given current trends, which 
can be divided into crop irrigated area proportion and irrigation norm. Subsequent to the beginning 
of 2010, food crop area ratio decreased sharply while the area of cash crops with high economic 
benefits has progressively increased (Table 6). Water-saving measures have developed slowly over 
time, and the primary pattern of irrigation has used traditional field canals. Advanced water-saving 
irrigation modes, such as micro-irrigation and greenhouse irrigation, only make up a smaller 
proportion. 

3.4.2 Scenario 2 


This setting was designed as an integrated agricultural water saving and water diversion scenario 
that enables a wider discussion of the integrated influence of the Dunhuang Oasis to the overall 
WCP. 

As agriculture is the largest water consumer within the Dunhuang Oasis, the implementation of 
water-saving measures in this context will be both indispensable and effective. Three concrete 
agricultural water saving-measures are currently utilized including (1) reducing irrigated area and 
planting area by 3600 hm? till 2020; (2) maintaining current food crop area, reducing the area of 
planted cotton and expanding the area of crops with high economic benefits (i.e., grapes and fruit 
trees); and (3) implementing water saving irrigation across the entire oasis, particularly by 
transforming the traditional irrigation mode to canal irrigation and micro-irrigation. 

Because of the water shortage caused by the agricultural water consumption, water saved by 
implementing water-saving measures will likely only maintain the water demand-supply balance 
under the current land use pattern and can hardly provide enough water for ecological services. In 
this context, the inter-basin water diversion project that encompasses the Sugan Lake Basin and the 
Danghe River Basin has been proposed to be developed over the next few years. This project plans 
to supply 0.835x10* m° of water into the Dunhuang Oasis from the Sugan Lake Basin. 

Analysis shows that the functions of agricultural and ecological subsystems will be modified 
under this scenario (Table 7) although remaining parameters will be the same as in Scenario 1. 
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Table 6 Parameters and Scenario 1 values 


Parameter 2010 2015 2020 2025 
Population growth rate (%bo) 5.10 5.10 5.10 5.10 
Urbanization rate (%) 50.00 59.00 65.00 65.00 
Rural residents living water consumption (m°/d) 65.00 65.00 65.00 65.00 
Urban residents living water consumption (m*/d) 130.00 130.00 130.00 130.00 
Large livestock water consumption (m*/d) 40.00 40.00 40.00 40.00 
Small livestock water consumption (md) 10.00 10.00 10.00 10.00 
Service industrial production growth rate (%) 7.12 7.12 11.38 9.96 
Service industrial water consumption (m*/(x10* CNY)) 14.00 11.30 9.30 9.30 
Industrial production growth rate (%) 7.50 8.70 9.40 10.50 
Industrial water consumption (m?/(x10* CNY)) 100.00 63.90 52.70 52.70 
irrigated area (x10* hm’) 3.52 3.52 3.52 3.52 
Food crop area proportion (%) 0.31 0.31 0.31 0.31 
Cotton crop area proportion (%) 53.02 46.57 46.57 46.57 
Grape area proportion (%) 19.13 21.13 21.13 21.13 
Melon area proportion (%) 2.63 2.63 2.63 2.63 
Fruit tree area proportion (%) 9.25 13.70 13.70 13.70 
Vegetable area proportion (%) 13.66 13.66 13.66 13.66 
Shelter forest area proportion (%) 2.00 2.00 2.00 2.00 
Food crop canal irrigation area proportion (%) 1.00 1.00 1.00 1.00 
Cotton micro-irrigation area proportion (%) 0.10 0.10 0.10 0.10 
Cotton canal irrigation area proportion (%) 0.90 0.90 0.90 0.90 
Vegetable canal irrigation area proportion (%) 0.91 0.91 0.91 0.91 
Vegetable micro-irrigation area proportion (%) 0.04 0.04 0.04 0.04 
Vegetable greenhouse micro-irrigation area proportion (%) 0.05 0.05 0.05 0.05 
Melon canal irrigation area proportion (%) 1.00 1.00 1.00 1.00 
Grape canal irrigation area proportion (%) 0.98 0.98 0.98 0.98 
Grape micro-irrigation area proportion (%) 0.02 0.02 0.02 0.02 
Fruit trees canal irrigation area proportion (%) 1.00 1.00 1.00 1.00 
Shelter forest canal irrigation area proportion (%) 1.00 1.00 1.00 1.00 
Food crop irrigation norm (m3/hm”) 6150.00 6150.00 6150.00 6150.00 
Cotton micro-irrigation norm (m*/hm?) 6750.00 6750.00 6750.00 6750.00 
Cotton canal irrigation norm (m*/hm?) 8250.00 8250.00 8250.00 8250.00 
Vegetable canal irrigation norm (m3/hm7) 10,200.00 10,200.00 10,200.00 10,200.00 
Vegetable micro-irrigation norm (m°/hm?) 7200.00 7200.00 7200.00 7200.00 
Vegetable greenhouse micro irrigation norm (m*/hm?) 6000.00 6000.00 6000.00 6000.00 
Melon irrigation norm (m?/hm?) 8250.00 7800.00 7800.00 7800.00 
Grape canal irrigation norm (m?/hm’) 9750.00 9750.00 9000.00 9000.00 
Grape micro-irrigation norm (m*/hm?) 6750.00 6750.00 6750.00 6750.00 
Shelter forest irrigation norm (m3/hm’) 6000.00 6000.00 6000.00 6000.00 
Water diversion (x10* m°) 0.00 0.00 0.00 0.00 
Flow into lower reaches (x10* m°) 4150.00 4150.00 4150.00 4150.00 
Crescent Spring water demand (x10* m°) 600.00 600.00 600.00 600.00 
Rural ecological water demand (x10* m°) 1677.00 1677.00 1677.00 1677.00 
Urban ecological water demand (x10* m°) 300.00 300.00 300.00 300.00 


3.4.3 Scenario 3 


This scenario comprises an ecological priority outcome that can be used to discuss how much 
irrigated land should be reclaimed to enable the ecological construction of the Dunhuang Oasis. 
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The national water resource strategy research report "Water Resource Configuration, Ecological 
Environment Construction, and Strategy of Sustainable Development Research in Northwest China 
(Water Volume)" (Chen et al., 2004) puts forward the ecological rehabilitation and its 
implementation across Northwest China. This report states that the water supply required for an 
ecosystem should not be less than 50% of the total regional water consumption. 

Based on Scenario 2, Scenario 3 uses the total water consumption of 2025 within the Dunhuang 
Oasis as the upper limit value in this ecological construction priority scenario. Thus, assuming that 
ecological water consumption is not less than 50% of the total water consumption, this analysis 
aims to determine how much irrigated land should be reduced, especially cash crops such as cotton, 
grapes and melon fields. This analysis does not take into account reductions in domestic and 
industrial water consumption because their proportions are so small. 


Table 7 Scenario 2 parameters 


Parameter 2010 2015 2020 2025 
Irrigated area (x10* hm?) 3.52 3.40 3.28 3.18 
Cotton crop area proportion (%) 53.02 43.07 39.20 37.20 
Grape area proportion (%) 19.13 24.60 27.50 29.50 
Fruit trees area proportion (%) 9.25 13.70 14.70 14.70 
Cotton micro irrigation area proportion (%) 0.10 0.29 0.39 0.40 
Cotton canal irrigation area proportion (%) 0.90 0.71 0.61 0.60 
Vegetable canal irrigation area proportion (%) 0.71 0.27 0.23 0.23 
Vegetable micro-irrigation area proportion (%) 0.22 0.54 0.54 0.54 
Vegetable greenhouse micro-irrigation area proportion (%) 0.07 0.19 0.23 0.23 
Grape canal irrigation area proportion (%) 0.80 0.43 0.40 0.40 
Grape micro irrigation area proportion (%) 0.20 0.57 0.60 0.60 
Water diversion (x10* m°) 0.00 0.00 8350.00 8350.00 
Flow into lower reaches (x10* m°) 4150.00 4150.00 8500.00 8500.00 
Groundwater recharge (x10* m°) 0.00 0.00 4000.00 4000.00 


4 Results and discussion 


4.1 Scenario 1 WCP 


The Scenario 1 prediction for 2025 implies that the total population of the Dunhuang Oasis will 
reach up to 19.39x 107; in this context, the urban population will increase from 9.40x10* in 2010 to 
12.59x10*4 in 2025 because of the urbanization. At the same time, domestic water consumption 
within the study area is likely to increase rapidly from a volume of 6.49 10° m? in 2010 to 7.54 10° 
m? in 2025. As a result of the booming tourism sector, the service industry is likely to consume 
almost 5.52* 10° m? of water in 2025, while industrial water consumption will increase sharply from 
10.88 10° m3 in 2010 to 22.60*10° m? in 2025. 


Table 8 Dunhuang Oasis water consumption pattern under Scenario 1 


Parameter 2010 2015 2020 2025 
Population (x10*) 17.97 18.44 18.91 19.40 
Urban population (x10*) 9.40 10.88 12.25 12.59 
Urban domestic water consumption (x10° m*) 4.46 5.16 5.81 5.97 
Rural domestic water consumption (x10° m°) 2.03 1.79 1.58 1.61 
Livestock water consumption (x10° m°) 2.25 2.46 2.68 2.92 
Industrial water consumption (x10° m°) 10.88 11.82 14.63 22.00 
Service industrial water consumption (x10° m°) 3.49 4.09 5.00 5.52 
Agricultural water consumption (x108 m°) 2.94 2.93 2.93 2.93 
Ecological water supply (x10’ m°) 2.57 2.57 2.57 2.57 


Total water consumption (x10* m°) 3.43 3.44 3.48 3.55 
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The volume of total available water within the Dunhuang Oasis is 4.02* 108 m3; in this context, 
the average surface water supply is 3.02x 108 m? while groundwater exploitation within the irrigated 
area remains at a maximum of 1.00x108 m?. This means that a volume of 0.415x108 m? must be 
removed from surface water to recharge the ecological system in the lower reaches, which would 
result in an ultimate usable surface water supply of 2.60x108 m°. The data presented in Table 
8summarize our SD simulation results. It is clear that total water consumption will increase from 
3.43108 m? in 2010 to 3.55x108 m? in 2025. These results suggest that if this oasis is to be 
maintained at its current developmental state, high groundwater exploitation intensity is essential 
even though the problem of groundwater depletion has already become increasingly serious. The 
world-famous natural landscape of the Crescent Spring, for example, has experienced a water level 
decline of 9-11 m over the last 40 years such that part of the lake bottom is now exposed during 
the low-water season (Zhou et al., 2007). We conclude that Scenario 1 is therefore inconsistent with 
sustainable development. 


4.2 Scenario 2 WCP 


The two most important ways to save water from agricultural irrigation are reducing the area treated 
and implementing conservation measurements. Thus, if a target reduction of 3600 hm? can be 
attained, then the volume of agricultural water that can be saved will be large. Agricultural water 
consumption will decrease by 0.27x108 m? given an average irrigation rate of 7500 m?/hm/. 
Subsequent to implementing water-saving measures and adjusting agricultural planting structures 
across the entire oasis, total agricultural water consumption is likely to decline significantly. The 
data summarized in Table 9 show that, compared to Scenario 1, the implementation of water-saving 


Table 9 Water consumption pattern comparisons 


Variable Year Scenario 1 Scenario 2 
2010 3.17 3.07 
Total water consumption! (x108 m°) 2015 oe 2.84 
2020 3.22 2001. 
2025 3.31 2.78 
2010 2.94 2.84 
2015 2.93 2.59 
Agricultural water consumption (x10* m°) 
2020 2.93 2.47 
2025 2.93 2.40 
2010 92.70 92.50 
; . 2015 92.00 91.10 
Agricultural water consumption (%) 
2020 90.70 89.20 
2025 89.00 86.30 
2010 3.43 3.32 
Total water consumption'+ecological water 2015 3.44 3.10 
supply (x10* m°) 2020 3.48 3.86 
2025 3.55 3.87 
2010 2.57 257 
2015 257: 2.57 
Ecological water supply (x107 m°) sae 2.57 10.92 
2025 2.57 10.92 
2010 7.49 7.74 
Ecological water supply (%) aie IA se 
2020 7.39 28.30 
2025 7.24 28.20 


Note: Total water consumption'=agricultural water subsystem+domestic water subsystem+industrial water subsystem+service water 
subsystem; agricultural water=agricultural water consumption/total water consumption!; ecological water=ecological water supply/(total 
water consumption!+ecological water supply). 
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measures within the Dunhuang Oasis are both very necessary and efficient as agricultural water 
consumption will decrease by almost 0.530108 m3. Agriculture remains the major consumer of 
water within the Dunhuang Oasis; thus, because of a large base, the proportion of total water 
consumption remains large. Models suggest that the proportion of water saved under Scenario 1 
will be 89.0% in 2025 while under the water saving Scenario 2, this value will be 86.3% (Table 9) 
(i.e., in order to emphasize water-saving effects, ecological water supply was subtracted from the 
denominator of total water consumption (Table 2)). The agricultural water that will be saved under 
this scenario can be assigned to domestic and industrial subsystems and could also contribute to 
reduce a negative groundwater balance. 

The additional water allocated from the inter-basin water diversion project also plays a very 
important role in alleviating the ecological water crisis within the oasis. Additional water diversion 
will increase ecological water supply in the oasis from 2.57x107 to 10.9210’ m? while the ratio of 
ecological water supply will also increase significantly. At the same time, according to the DBP, 
groundwater exploitation will decrease from the current 10.0x107 to 6.80107 m? in the future. The 
increase of ecological water can recharge the groundwater level within the Dunhuang Oasis and 
improve the state of the lower reaches of West Lake Wetland. 


4.3 Scenario 3 WCP 


The ecological water supply under Scenario 2 is likely to reach a maximum volume of 10.9x107m? 
in 2025 while total water consumption is 3.87x108 m?. The assumption above states that the 
ecological water supply should not be less than 50% of total water consumption while the actual 
supply in this case should reach a goal of 1.94x108 m? such that agricultural water consumption 
will be decreased by 0.84% 108 m°. The results of our simulations show that if the overall proportions 
of crop structure remain unchanged, then total irrigated area will be decreased to 2.06 x104 hm? by 
2025 (Fig. 10). Areas of cotton, grapes, fruit trees and vegetables have also greatly decreased and 
play a critical role in Scenario 3 (Fig. 8). 


1.80 
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1.60 = Melons m Tree-grasses ™ Food crops 


1.40 
1.20 
1.00 
0.80 


Area (x 104 hm’) 


0.60 
0.40 
0.20 


0.00 


Scenario 1 Scenario 2 Scenario 3 


Fig. 10 Crop areas under the three scenarios 


5 Conclusions 


We analyzed WCP within the Dunhuang Oasis under three scenarios using a system dynamics 
model. The results of this study reveal that agriculture is the primary consumer of water within the 
Dunhuang Oasis (>90%) while domestic, industrial and service industrial water consumption all 
together play a small role in Dunhuang Oasis WCP. As the total water consumption of these three 
accounted for 6.7%, 7.4%, 8.5% and 9.7% in 2010, 2015, 2020 and 2025, respectively, limited 
savings are likely to be found from these sections. Data show that under Scenario 1, the present- 
day baseline, total water consumption within the Dunhuang Oasis is likely to progressively increase 
as agricultural uses remain quite high. Depletion at this level will lead to worsening groundwater 
losses and will threaten the ecological security of the oasis. Modeling results suggest that 
implementing water-saving measures will lead to savings of almost 5.30x107 m3. A substantial 
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agricultural base will lead to large total water consumption; thus, by 2025, the proportion of water 
consumed by agriculture will decrease by only 2.7% from 89.0% under Scenario 1 to 86.3% under 
Scenario 2. The implementation of water-saving measures cannot be ignored; agricultural water 
saved in this way can be re-assigned to the domestic and industrial subsystems and can also 
contribute to mitigating groundwater depletion. The additional water supply available from the 
inter-basin diversion project will increase ecological water demand from 2.57* 107 to 10.9210? m? 
while the water supply in this regard will also increase significantly. Additional water allocated 
from the inter-basin water diversion project will also play a very important role in alleviating the 
oasis ecological water crisis. Thus, if regional ecological construction is regarded as the most 
important development strategy within the Dunhuang Oasis, total irrigated area should be decreased 
to 2.06104 hm? by 2025. 
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